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ADDITIONAL ASPECTS OF ELASTOHYDRODYNAMIC LUBRICATION 


by Bernard J. Hamrock 
Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 

ABSTRACT 

Professor Cheng is to be congratulated for doing an outstanding Job 
in presenting an up-to-date review of the varying aspects of elastohydro- 
dynamic lubrication. In the present paper I would merely like to add 
some material which was not covered in Professor Cheng's paper. In par- 
ticular 1 would like to present some recent work on elastohydrodynamic 
lubrication of materials of low elastic modulus as well as work on hydro- 
dynamic lubrication. Both these topics will be applicable for contacts 
with any ellipticity parameter (ranging from a circular contact to a line 
contact). These two studies will be combined with some previous work to 
supply film thickness contours for the different regimes of lubrication 
as a function of the elasticity and viscosity parameters. From these 
studies the lubrication regime and the minimum film thickness within the 
contact can be easily determined. 
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ADDITIONAL ASPECTS OF ELASTOHYDRODYNAMIC LUBRICATION 


by Bernard J. Hamrock 
Lewis Research Center 
INTRODUCTION 

Figure l points out the difference between hydrodynamic and elasto- 
hydrodynamic lubrication (EHL). Hydrodynamic lubrication normally deals 
with con formal contacts while elastohydrodynamic lubrication deals with 
noncon formal contacts. Elastohydrodynamic lubrication can be divided 
into materials of low and high elastic modulus. Note the maximum pres- 
sure is three orders of magnitude less for materials of low elastic mod- 
ulus as compared to that ot high elastic modulus. Also, the minimum film 
thickness is the smallest tor EHL ot materials of high elastic modulus 
and next smallest for EHL ot materials of low elastic modulus and the 
largest tor hydrodynamic lubrication contacts. 

Lubricant in elastohydrodynamic lubrication of high elastic modulus 

materials enters a typical contact with atmospheric pressure, is subjected 

2 

to a maximum pressure of approximate Iv J00 01)0 lb/in. “, traverses a passage 
one thousand times longer than its height, and Is then ejected into the 
atmosphere within a time* span ot about 10 - * or 10"^ second. The viscous 
character of the fluid while passing through the contact changes drasti- 
cally, going from an easy flowing liquid to a pseudosolid and back to an 
easy flowing liquid In less than a fraction of a millisecond. 

Lubricant in elastohydrodynamic lubrication of low elastic modulus 
materials enters a typical contact with atmospher ic pressure and is sub- 
jected only to a maximum pressure of J^IO* - psi due to the elastomeric 
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nature ot tlie material. For such materials the distortions are large 
even with light loads. Another feature of EHL of low elastic materials 
is the negligible efiect of pressure on the viscosity of the lubricating 
fluid. Therefore the viscous character of the fluid does not change 
while passing through the contact. 

NOMENCLATURE 

a semimajor axis of contact ellipse 

b semi mi nor axis of contact ellipse 

E modulus of elasticity 


> 



elliptical integral of second kind 
F normal applied load 

G dimensionless material parameter, aE' 

H dimensionless film thickness, h/R 

x 

H dimensionless film thickness, H(W/U)“ 
h film thickness 

k elliptlcity parameter, a/b 

m dimensionless inlet distance 

m* dimensionless inlet distance at boundary between fully flooded and 
starved conditions 
R effective radius 

r radius of curvature 

2 

U dimensionless speed parameter, un^/E'R^ 

u surface velocity in rolling direction, (u. + u_)/2 

A B 
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W dimensionless load parameter, F/E'R^ 

a pressure-viscosity constant 

B R/R 

y x 

Hq atmospheric viscosity 
v Poisson's ratio 



Subscripts : 

A solid A 

B solid B 

E elastic 

F fully flooded conjunction 

1 isoviscous 

min minimum 

R rigid 

S starved conjunction 

V viscous 

x,y coordinate system 

EHL FOR MATERIALS OF LOW ELASTIC MODULUS 
The work, presented by Professor Cheng related only to materials of 
high elastic modulus (e.g., steel). The work that I would like to present 
in this section is for materials of low elastic modulus (e.g., nitrile 
rubber). Engineering applications in which elastohydrodynamic lubrication 
is important for low- elastic-modulus materials include seals, human joints, 
tires, and elastomeric-material machine elements. 
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The problem of fully flooded line contacts has been solved theoret- 
ically for low -el as tic -mate rials by Herrebrugh (1), Dowson and Swales (2), 
and Boglin and Archard (3). The solutions of references 1 and 2 were ob- 
tained numerically and are based on simultaneous solutions of the hydro- 
dynamic and elasticity equations; the analytical solution of reference 3 
relied on the assumption of a simplified form for the film shape in the 
contact region. Biswas and Snidle (4) used the approach of reference 3 
to solve the point-contact situation. Reference 5 presents, to the best 
of the author's knowledge, the first complete numerical solution of the 
problem of fully flooded isothermal elastohydrodynamic lubrication of 
elliptical contacts for low-elastic-modulus materials. Reference 6 ex- 
tends the work of reference 5 by studying the effect of lubricant starva- 
tion on pressure and film thickness within the conjunction. This section 
of the present paper will utilize the work of references 5 and f >. 

Reference 3 produces the following dimensionless minimum film thick- 
ness for fully flooded elliptical contacts for materials of low elastic 
modulus 


(H . 

nun ,F 


- 7.43(1 - 0.85 e -0,31 k )U°* 65 W _0,21 


( 1 ) 


l.E 


The subscript 1,E refers to the isoviscous elastic lubrication condi- 
tion which is the condition that exists in dealing with materials of low 
elastic modulus. 

It is interesting to compare the equation for materials of low elastic 
modulus (eq. (1)) with the corresponding equation found in reference 7 for 
materials of high elastic modulus, namely: 
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H irin,F 


, -0.68 k. l .0.b8 tJ -0. 07 3^0.49 

* 3.b3(l - e )U W G 


V.E 


( 2 ) 


File powers of U in equations (1) and (2) are quite similar, but the 
power of W is much more significant for low -e last lc -modulus materials. 
The expressions showing the effect of the ellipticity parameter is of ex- 
ponential form in both equations, but with quite different constants. 

A major difference in equation (1) and (2) is the absence of a mate- 
rial parameter in the expression lor the minimum film thickness tor low- 
elast ic-modulus materials. The reason for this is the negligible effect 
of pressure on the viscosity of the lubricating fluid for low-elastic- 
modulus materials. 

When the elliptical contact becomes starved, the dimensionless mini- 
num film thickness for materials of low elastic modulus (6) can be written 
as 


(H . c ) 

min ,S 


1 »E 


(H , — \ 

mln * F 1 ,E \m* - I' 


0.22 


(3) 


where m is the distance of the inlet miniscus from the center of the 
contact, and m* is the inlet distance required tor achieving tiie fully 
flooded condition. From reference b n* can be expressed as 


m 


* . 



0.1b 


{H c } 
min ,1 


I.E 


( 4 ) 


where b is the semiminor axis of the elliptical contact. 

To make it easier to calculate the semiminor axis of the contact 
ellipse (b), the elliptical integral of the second kind (,* T ) , and the 
ellipticity parameter, the approximate expressions developed in refer- 
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ence 8 will be used and are given below as 


bis 


irkE 




1/3 


(5) 


3R 

^ * 1 + 5T 
y 


/R X 0 - 64 

■ te) 


( 6 ) 


(7) 


where 


r Ax r Bx 


-L.+ JL 

r Ay r By 


( 8 ) 


(9) 


The approximate expressions (eqs. (5) to (7)) enable one to easily calcu- 
late these terms within 3 percent accuracy without resorting to charts or 
numerical methods. 

To explain more fully what happens to the film thickness in going 
from a fully flooded to a lubricant starvation condition for materials of 
low elastic modulus, figures 2 and 3 are presented. Figure 2 represents 
a fully flooded condition and figure 3 a severely starved condition. In 
these figures the symbol + indicates the center of the Hertzian contact. 
Because of the way the coordinates are made dimensionless, the actual con- 
tact ellipse becomes a Hertzian circle regardless of the ellipticity param- 
eter. The Hertzian contact circle is shown in these figures by asterisks. 
At the top of these figures the contour label and its corresponding value 
is given. The inlet region is to the left and the exit region is to the 


right. 
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It is observed in these figures that the central portion of the 
film-thickness contours has become parallel for the severely starved con- 
dition (fig. 3). Hie minimum film thickness area is closer to the exit 
region for the severely starred condition. Note also that the values of 
the film thickness contours for the severely starved condition (fig. 3) 
are much lower than those of the fully flooded condition (fig. 2). 

DIMENSIONLESS GROUPING 


The dimensionless group {H, U, W, G , and k} used in equations (1) 
and (2) have served a very useful rpose in aiding the understanding of 
the results found in references 5 to 7. Several authors (e.g., refs. 9 
and 10) have noted that the examination of the dimensionless representa- 
tion of governing equations shows the above set of dimensionless group 
cun be reduced by one without any loss of generality. Although Johnson's 
paper (10) does not consider elliptical contacts, it does state what the 
nondimensional parameters would be, namely: 

Dimensionless film thickness 


min 


H mln (u)’ 


Dimensionless viscosity parameter 


GW' 

l i " :r 


Dimensionless elasticity parameter 


U 


8/3 


*3 2 

U 


( 10 ) 


( 11 ) 


( 12 ) 


The ellipticity parameter (k) remains a parameter as discussed in equa- 
tion (7). Therefore, the reduced dimensionless group is {H , g , g 

min 1 J 


and k }. 
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FILM THICKNESS IN ELLIPTICAL CONTACTS FOR FOUR LUBRICATION REGIMES 
Reference 11 makes use of the reduced dimensionless grouping to give 
Che film thickness equations for four fluid film lubrication regimes 
found in elliptical contacts. These regimes are distinguishable by the 
influence or lack of influence of elastic and viscous effects. The film 
thickness equations for the respective regimes come from earlier theoret- 
ical studies on elastohydrodynamic (5 and 7) and hydrodynamic (12) lubri- 
cation of conjunctions ot elliptical form. 

ISOVISCOUs RIGID REGIME 


The influence of geometry on the hydrodynamic film separating two 
rigid solids was investigated in reference 12. It was found that the 
minimum film thickness had the speed, viscosity, and load dependence as 
the classical Kapitza (13) solution. However, the incorporation of 
Reynolds condition resulted in an additional geometry effect. Therefore, 
from reference 12 the dimensionless minimum film thickness for the iso- 
viscous rigid lubrication regime can be written as 

{H , ) - 128dS'|i.o« tan' 1 f + 0.13 (13) 

" in I.R l - j 

whe re 


6 



_k\l /0 - 6 * 

1.03/ 


(14) 



(15) 


In equation (13) note that the dimensionless film thickness is strictly a 
function of the geometry of the contact (R^/R^). 
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VISCOUS RIGID REGIME 


Making use of Blok (14) while adding a geometry effect for an ellip- 
tical contact as discussed in reference 11, the minimum film thickness 
for the viscous rigid regime can be written as 


* 2/3 

1H : ) * 1-bb gT (1 ' 

mi" VK 1 


—0 . 68 k . 
e ) 


(lb) 


ISOVISCOUS ELASTIC REGIME 


Expressing equation (1) in terms of the reduced dimensionless group 

(H , . g, . and k) 1 rom reference 11 the dimensionless minimum film thick- 

min* & 1 ’ 

ness for the isoviscous elastic regime can be written as 


(H . ) « 8.70 g°’ b7 (l - 0.85 e -0 * 31 k ) 

min 3 

1 9 t. 


(17) 


VISCOUS ELASTIC 


Expressing equation (2) in terms of the reduced dimensionless group 

(H . . g, . g,. and 0 the minimum film thickness for the viscous elastic 
min’ “l* °3* 

regime can be written as 


, ... 0.49 0.17,. -0.68 k, 

‘"min' * 3 ' 45 *1 *3 (l ' e > 

V 


(18) 


Having defined the dimensionless minimum film thickness for the four 
lubrication regimes in equations (13), (lb), (17), and (18) these equa- 
tions were used to develop the dimensionless film thickness contours for 
the different regimes of lubrication. Reference 11 shows these maps on a 
log-log grid of dimensionless viscosity parameter and dimensionless elas- 
ticity parameter for ellipticity parameters of 1, 2, 3, 4, and 6. Fig- 
ure 4 shows a sample of tnese results for an ellipticity parameter of 1. 
Tiie four lubrication regimes are clearly indicated in this figure. In 
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this figure we >ee that the transition between the viscous rigid and vis- 
cous elastic is nearly smooth whereas the transition from the isoviscous 
elastic and viscous elastic is not as good. One right speculate that the 
transition t rom the isoviscous elastic and viscous elastic i..ignt produce 
sligntly lower dimensionless film thicknesses than those in figure 4. 
Similar results are shown for other ellipticity parameters in refer- 
ence 11. 

Theretore, given an ellipticity parameter (k) , the value of the vis- 
cosity parameter and the elasticity parameter (gj) , one can deter- 

mine tne lubricating regime as well as an approximate value of the dimen- 
sionless minimum film thickness. knowing the lubrication regime a more 

accurate value ot H . can be obtained using the appropriate dimension- 

min 

less minimum tilm thickness equation (one of tire following equations: 
(13). (lb), (17), or (18)). 
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ELASTIC EFFECTS PREDOMINATE ELASTIC AND VISCOUS EFFECTS 

BOTH IMPORTANT 

Figure 1. - Hydrodynamic and elastohydrodynamic lubrication 





Figure 3. - Contour plot of dimensionless film thickness is a sewrly 
starved conjunction for materials of low elastic modulus. 



